A Novel Subcluster of Closely Related Bacillus Phages with Distinct Tail Fiber/Lysin Gene Combinations.

Bacteriophages (phages) are the most numerous entities on Earth, but we have only scratched the surface of describing phage diversity. We isolated seven Bacillus subtilis phages from desert soil in the southwest United States and then sequenced and characterized their genomes. Comparative analyses revealed high nucleotide and amino acid similarity between these seven phages, which constitute a novel subcluster. Interestingly, the tail fiber and lysin genes of these phages seem to come from different origins and carry out slightly different functions. These genes were likely acquired by this subcluster of phages via horizontal gene transfer. In conjunction with host range assays, our data suggest that these phages are adapting to hosts with different cell walls.

Genomic and phenotypic signatures of bacteriophage coevolution with the phytopathogen Pseudomonas syringae.

The rate and trajectory of evolution in an obligate parasite is critically dependent on those of its host(s). Adaptation to a genetically homogeneous host population should theoretically result in specialization, while adaptation to an evolving host population (i.e., coevolution) can result in various outcomes including diversification, range expansion, and/or local adaptation. For viruses of bacteria (bacteriophages, or phages), our understanding of how evolutionary history of the bacterial host(s) impacts viral genotypic and phenotypic evolution is currently limited. In this study, we used whole genome sequencing and two different metrics of phage impacts to compare the genotypes and phenotypes of lytic phages that had either coevolved with or were repeatedly passaged on an unchanging (ancestral) strain of the phytopathogen Pseudomonas syringae. Genomes of coevolved phages had more mutations than those of phages passaged on a constant host, and most mutations were in genes encoding phage tail-associated proteins. Phages from both passaging treatments shared some phenotypic outcomes, including range expansion and divergence across replicate populations, but coevolved phages were more efficient at reducing population growth (particularly of sympatric coevolved hosts). Genotypic similarity correlated with infectivity profile similarity in coevolved phages, but not in phages passaged on the ancestral host. Overall, while adaptation to either host type (coevolving or ancestral) led to divergence in phage tail proteins and infectivity patterns, coevolution led to more rapid molecular changes that increased bacterial killing efficiency and had more predictable effects on infectivity range. Together, these results underscore the important role of hosts in driving viral evolution and in shaping the genotype-phenotype relationship.

Forty Years without Family: Three Novel Bacteriophages with High Similarity to SPP1 Reveal Decades of Evolutionary Stasis since the Isolation of Their Famous Relative.

SPP1, an extensively studied bacteriophage of the Gram-positive Bacillus subtilis, is a model system for the study of phage-host interactions. Despite progress in the isolation and characterization of Bacillus phages, no previously fully sequenced phages have shared more than passing genetic similarity to SPP1. Here, we describe three virulent phages very similar to SPP1; SPP1 has greater than 80% nucleotide sequence identity and shares more that 85% of its protein coding genes with these phages. This is remarkable, given more than 40 years between the isolation of SPP1 and these phages. All three phages have somewhat larger genomes and more genes than SPP1. We identified a new putative gene in SPP1 based on a conserved sequence found in all phages. Gene conservation connotes purifying selection and is observed in structural genes and genes involved in DNA metabolism, but also in genes of unknown function, suggesting an important role in phage survival independent of the environment. Patterns of divergence point to genes or gene domains likely involved in adaptation to diverse hosts or different environments. Ultimately, comparative genomics of related phages provides insight into the long-term selective pressures that affect phage-bacteria interactions and alter phage genome content.

Complete Genome Sequence of Finnry, a Subcluster L3 Mycobacteriophage from Charleston, South Carolina.

Subcluster L3 bacteriophage Finnry was isolated from soil collected in Charleston, South Carolina, using Mycobacterium smegmatis mc2155 as a host. The genome of this temperate siphovirus is 75,632 bp long (130 predicted protein-coding genes, 9 tRNAs, and no transfer-messenger RNAs), and BLASTn alignment revealed 99.86% identity with the genome of L3 mycobacteriophage Samty.

Complete Genome Sequences of Two Temperate Bacillus subtilis Phages Isolated at Tumamoc Hill Desert Laboratory.

Bacteriophages are important in structuring bacterial communities, including desert soils dominated by Bacillus species. Here, we describe two genetically similar temperate phages isolated on a Bacillus subtilis strain from soil in Tucson, Arizona. Their double-stranded DNA (dsDNA) genomes contain 98 and 102 genes, with a set of 4 genes being found in only one phage.

Isolation and Annotation of the Genome Sequences of Bacteriophages InvictusManeo (Subcluster K5) and Netyap (Subcluster L2).

The mycobacteriophages InvictusManeo (K5 subcluster) and Netyap (L2 subcluster) were isolated from soils in Cullowhee Creek, Cullowhee, North Carolina. Both exhibit Siphoviridae morphology and infect Mycobacterium smegmatis mc2155. The InvictusManeo genome is 61,147 bp and contains 96 predicted protein-coding genes, whereas the Netyap genome is 76,366 bp with 131 predicted protein-coding genes.

Comparative Genomics of Six Lytic Bacillus subtilis Phages from the Southwest United States.

Background: Despite their importance to microbial dynamics involving Bacillus subtilis, we have a limited understanding of the diversity of phages that can lyse this model organism. Materials and Methods: Phages were isolated from soil samples collected from various sites in the southwest U.S. deserts on a wild B. subtilis strain. Their genomes were assembled, characterized, and bioinformatically compared. Results: Six Siphoviruses with high nucleotide and amino acid similarity to each other (>80%) but very limited similarity to phages currently in GenBank were isolated. These phages have double-stranded DNA genomes (55,312 to 56,127 bp) with 86-91 putative protein coding genes, and a low GC content. Comparative genomics reveal differences in loci encoding proteins that are putatively involved in bacterial adsorption with evidence for genomic mosaicism and a possible role for small genes. Conclusions: A comparative approach provides insights into phage evolution, including the role of indels in protein folding.

Complete Genome Sequences of Four Phages of the Horse Chestnut Phyllosphere.

Bacteriophages play important roles in determining bacterial communities, including plant microbiota. Here, we describe four lytic phages, three Siphoviridae and one Podoviridae, isolated from four different bacterial species found on the leaves of horse chestnut trees. Their double-stranded DNA (dsDNA) genomes range from 39,095 to 46,062 bp and contain 51 to 70 genes.

Genome Sequences of Erwinia Phyllophages AH04 and AH06.

Although crucial in shaping bacterial communities, few bacteriophages of the phyllosphere have been described. We provide genome data for two Myoviridae phages, AH04 and AH06, isolated on Erwinia billingiae strains. AH04 shares limited genetic similarity with previously described phages, while AH06 shares over 75% similarity with various Erwinia phages.

Complete Genome Sequence of the Pantoea Phage AH07.

Bacteriophages of the phyllosphere have not been extensively described, despite their role in bacterial communities on this plant organ. Here, we describe a temperate Pantoea phage, AH07, that was isolated from the leaves of horse chestnut trees. The 37,859-bp linear double-stranded DNA genome contains 58 putative genes, including an integration cassette.

Complete Genome Sequence of the Cluster B4 Mycobacteriophage Lolalove, Isolated in Charleston, South Carolina.

Lolalove, a B4 subcluster soil bacteriophage of Mycobacterium smegmatis, was isolated in Charleston, SC. It possesses a 71,111-bp linear double-stranded DNA (dsDNA) genome with 99 protein-coding genes and a GC content of 68.9%. genome BLASTn alignments indicate high sequence identity to the related B4 subcluster M. smegmatis phages BrownCNA, Mithril, and Hangman.

Genome Sequences of the Mycobacteriophages Awesomesauce and LastJedi.

Bacteriophages Awesomesauce and LastJedi infect Mycobacterium smegmatis mc2155. While the Awesomesauce genome is 57,054 bp with 94 protein-coding genes, the LastJedi genome is 55,149 bp with 94 protein-coding genes. Nucleotide sequence comparison in Phamerator detected synteny between Awesomesauce gp49 to gp61 and singleton LilSpotty. Whole-genome BLASTn alignments revealed that LastJedi strongly resembles Clifton (99.41% identity).

Complete Genome Sequence of the Cluster F1 Mycobacteriophage KingMidas.

Subcluster F1 bacteriophage KingMidas was isolated from soil collected in Providence, Rhode Island, using Mycobacterium smegmatis mc2155 as the host. The genome is 57,386 bp and contains 105 predicted protein-coding genes but no transfer-messenger RNAs or tRNAs. This siphovirus has an icosahedral head, with a genome 99.1% identical to that of F1 mycobacteriophage Scottish.

Genome Sequences of Arthrobacter sp. Phages Inspire2, Ronnie, Hunnie, DeWayne, CGermain, Copper, Azathoth, and Arby.

Eight siphoviral phages isolated from various soil types and locations in southwestern Pennsylvania using Arthrobacter sp. strain ATCC 21022 were sequenced. The phages all have relatively small genomes, with each genome containing 15,556 bp. All 8 phages are closely related to previously described cluster AN Arthrobacter phages (K. K. Klyczek, J. A. Bonilla, D. Jacobs-Sera, T. L. Adair, et al., PLoS One 12:e0180517, 2017, https://doi.org/10.1371/journal.pone.0180517; J. Y. Lee-Soety, S. Bhatt, T. L. Adair, J. A. Bonilla, et al., Genome Announc 5:e01092-17, 2017, https://doi.org/10.1128/genomeA.01092-17).

Multi-Institutional, Multidisciplinary Study of the Impact of Course-Based Research Experiences.

Numerous national reports have called for reforming laboratory courses so that all students experience the research process. In response, many course-based research experiences (CREs) have been developed and implemented. Research on the impact of these CREs suggests that student benefits can be similar to those of traditional apprentice-model research experiences. However, most assessments of CREs have been in individual courses at individual institutions or across institutions using the same CRE model. Furthermore, which structures and components of CREs result in the greatest student gains is unknown. We explored the impact of different CRE models in different contexts on student self-reported gains in understanding, skills, and professional development using the Classroom Undergraduate Research Experience (CURE) survey. Our analysis included 49 courses developed and taught at seven diverse institutions. Overall, students reported greater gains for all benefits when compared with the reported national means for the Survey of Undergraduate Research Experiences (SURE). Two aspects of these CREs were associated with greater student gains: 1) CREs that were the focus of the entire course or that more fully integrated modules within a traditional laboratory and 2) CREs that had a higher degree of student input and results that were unknown to both students and faculty.

Testing hypotheses for the presence of tRNA genes in mycobacteriophage genomes.

The presence of tRNA genes in bacteriophages has been explained on the basis of codon usage (tRNA genes are retained in the phage genome if they correspond to codons more common in the phage than in its host) or amino acid usage (independent of codon, the amino acid corresponding to the retained tRNA gene is more common in the phage genome than in the bacterial host). The existence of a large database of sequenced mycobacteriophages, isolated on the common host Mycobacterium smegmatis, allows us to test the above hypotheses as well as explore other hypotheses for the presence of tRNA genes. Our analyses suggest that amino acid rather than codon usage better explains the presence of tRNA genes in mycobacteriophages. However, closely related phages that differ in the presence of tRNA genes in their genomes are capable of lysing the common bacterial host and do not differ in codon or amino acid usage. This suggests that the benefits of having tRNA genes may be associated with either growth in the host or the ability to infect more hosts (i.e., host range) rather than simply infecting a particular host.

Complete Genome Sequences of Five Bacteriophages That Infect Rhodobacter capsulatus.

Five bacteriophages that infect the Rhodobacter capsulatus strain YW1 were isolated from stream water near Bloomington, Illinois, USA. Two distinct genome types are represented in the newly isolated bacteriophages. These genomes are different from other bacteriophage genomes previously described.

Cluster M mycobacteriophages Bongo, PegLeg, and Rey with unusually large repertoires of tRNA isotypes.

UNLABELLED: Genomic analysis of a large set of phages infecting the common host Mycobacterium smegmatis mc(2)155 shows that they span considerable genetic diversity. There are more than 20 distinct types that lack nucleotide similarity with each other, and there is considerable diversity within most of the groups. Three newly isolated temperate mycobacteriophages, Bongo, PegLeg, and Rey, constitute a new group (cluster M), with the closely related phages Bongo and PegLeg forming subcluster M1 and the more distantly related Rey forming subcluster M2. The cluster M mycobacteriophages have siphoviral morphologies with unusually long tails, are homoimmune, and have larger than average genomes (80.2 to 83.7 kbp). They exhibit a variety of features not previously described in other mycobacteriophages, including noncanonical genome architectures and several unusual sets of conserved repeated sequences suggesting novel regulatory systems for both transcription and translation. In addition to containing transfer-messenger RNA and RtcB-like RNA ligase genes, their genomes encode 21 to 24 tRNA genes encompassing complete or nearly complete sets of isotypes. We predict that these tRNAs are used in late lytic growth, likely compensating for the degradation or inadequacy of host tRNAs. They may represent a complete set of tRNAs necessary for late lytic growth, especially when taken together with the apparent lack of codons in the same late genes that correspond to tRNAs that the genomes of the phages do not obviously encode. IMPORTANCE: The bacteriophage population is vast, dynamic, and old and plays a central role in bacterial pathogenicity. We know surprisingly little about the genetic diversity of the phage population, although metagenomic and phage genome sequencing indicates that it is great. Probing the depth of genetic diversity of phages of a common host, Mycobacterium smegmatis, provides a higher resolution of the phage population and how it has evolved. Three new phages constituting a new cluster M further expand the diversity of the mycobacteriophages and introduce novel features. As such, they provide insights into phage genome architecture, virion structure, and gene regulation at the transcriptional and translational levels.

Stability of pollen-ovule ratios in pollinator-dependent versus autogamous Clarkia sister taxa: testing evolutionary predictions.

It has been proposed that natural selection should favor distinct temporal patterns of sex allocation in selfing vs pollinator-dependent taxa. In autogamous selfers in which pollen receipt is highly reliable, selection should favor genotypes that maintain low and stable pollen to ovule (P : O) ratios throughout flowering. By contrast, in outcrossers the optimum P : O ratio of an individual's flowers will depend on pollinator abundances and mating opportunities, both of which may vary over time. In this case, selection may favor temporal variation among flowers in the P : O ratio. An opposing prediction is that selfing taxa will be developmentally more unstable than outcrossers because of lower homeostasis caused by high homozygosity. We compared temporal changes in the P : O ratio in two pairs of sister taxa in the genus Clarkia. We examined hundreds of glasshouse-raised maternal families representing three wild populations each of the outcrossing, insect-pollinated Clarkia unguiculata, the facultatively autogamous Clarkia exilis and the outcrossing and selfing subspecies of Clarkia xantiana: ssp. xantiana and parviflora, respectively. Temporal change in the P : O ratio was significantly greater in both outcrossers than in their selfing sister taxa, although the proportional changes in the P : O ratio (relative to the first bud produced) did not differ significantly between sister taxa (0.07 < P < 0.10). Our results provide partial support for the hypothesis that the P : O ratio is more stable in selfing than in outcrossing taxa and reject the hypothesis that selfers are less stable.

Size-dependent pollen:ovule ratios and the allometry of floral sex allocation in Clarkia (Onagraceae) taxa with contrasting mating systems.

Multiple field populations of two pairs of diploid sister taxa with contrasting mating systems in the genus Clarkia (Onagraceae) were surveyed to test predictions concerning the effects of resource status, estimated as plant size, on pollen and ovule production and on the pollen:ovule (P:O) ratio of flowers. Most theoretical models of size-dependent sex allocation predict that, in outcrossing populations, larger plants should allocate more resources to female function. Lower P:O ratios in larger plants compared to smaller plants have been interpreted as supporting this prediction. In contrast, we predicted that P:O ratio should not vary with plant size in predominantly selfing plants, in which each flower contributes to reproductive success equally through male and female function. We found that, in all four taxa, both ovule and pollen production per flower usually increased significantly with plant size and that the shape of this relationship was decelerating. However, ovule production either decelerated more rapidly than or at the same rate as pollen production with plant size. Consequently,the P:O ratio increased or had no relationship with plant size. This relationship was population-specific (not taxon-specific) and independent of the mating system. Possible explanations for the increasing maleness with plant size are discussed.